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ABSTRACT 



In order to explore the properties of extreme outer stellar disks, we obtained 
ultra-deep V and GALEX UV images of 4 dwarf irregular galaxies and one Blue 
Compact Dwarf galaxy and ultra-deep B images of 3 of these. Our V-band sur- 
face photometry extends to 29.5 magnitudes arcsec"^. We convert the FUV and 
V-band photometry, along with Ha photometry obtained in a larger survey, into 
radial star formation rate profiles that are sensitive to timescales from 10 Myrs to 
the lifetime of the galaxy. We also obtained H I -line emission data and compare 
the stellar distributions, surface brightness profiles, and star formation rate pro- 
files to Hi -line emission maps, gas surface density profiles, and gas kinematics. 
Our data lead us to two general observations: First, the exponential disks in these 
irregular galaxies are extraordinarily regular. We observe that the stellar disks 
continue to decline exponentially as far as our measurements extend. In spite of 
lumpiness in the distribution of young stars and H I distributions and kinematics 
that have significant unordered motions, sporadic processes that have built the 
disks — star formation, radial movement of stars, and perhaps even perturbations 
from the outside — have, nevertheless, conspired to produce standard disk profiles. 
Second, there is a remarkable continuity of star formation throughout these disks 
over time. In four out of five of our galaxies the star formation rate in the outer 
disk measured from the FUV tracks that determined from the l^-band, to within 
factors of 5, requiring star formation at a fairly steady rate over the galaxy's 
lifetime. Yet, the Hi surface density profiles generally decline with radius more 
shallow ly than the stellar light, and the gas is marginally gravitationally stable 
against collapse into clouds. Outer stellar disks are challenging our concepts of 
star formation and disk growth and provide a critical environment in which to 
understand processes that mold galaxy disks. 

Subject headings: galaxies: irregular — galaxies: star formation — galaxies: 
individual (DDO 53, DDO 86, DDO 133, NGC 4163, IZw 115) — galaxies: 
structure 



1. Introduction 

Outer parts of galaxy disks present an extreme environment for star formation. Dwarf 
irregular (dim) galaxies already challenge models of star formation because of their low 
gas densities even in the central regions (Hunter & Plummer 1996; Meurer et al. 1996; 
van Zee et al. 1997; Hunter, Elmegreen, & Baker 1998). Outer parts of dim galaxies with 
their exceptionally low gas densities, therefore, present a particularly difficult test of our 
understanding of the cloud/star formation process. Furthermore, knowledge of the star 
formation histories of outer stellar disks provides information on the way in which the disk 
has grown with time and tests current models for disk growth (see, for example, Azzollini et 
al. 2008). 

Stars have formed even to very low disk surface brightness levels. For example, Saha 
et al. (2010) have detected an exponential stellar disk in the LMC that extends to 12 disk 
scale lengths, the equivalent of an J-band surface brightness of 34 mag arcsec"^. Although 
the number of young stars relative to older stars declines with radius, stars as young as a 
few Gyr are found to about 8.3 disk scale lengths. In the late-type spiral M33, Barker et al. 
(2011) have found stars younger than 0.5 Gyr at 4 and 5 disk scale lengths, although the bulk 
of the stars are older. And in M31, Richardson et al. (2008) find evidence for star formation 
over the past 0.25-2 Gyr as far out as 56 kpc. Furthermore, satellite ultraviolet imaging 
has revealed extended ultraviolet disks in many spiral galaxies that indicate on-going star 
formation in outer disks (Thilker et al. 2007). For example, GALEX images of M83 revealed 
UV-bright stellar complexes to 7 kpc, about 4 times the radius at which most of the H ll 
regions were detected (Thilker et al. 2005). 

Given the importance of outer stellar disks to our understanding of galaxy formation 
and evolution, we have conducted a program to detect and measure the extreme outer stellar 
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disks of a sample of dim galaxies through ultra-deep imaging in the optical and ultraviolet. 
In addition we have mapped the galaxies in H I -line emission. We use these data to examine 
the structure, stellar populations, and star formation activity in extreme outer dim disks. 

Our program objects were selected from a large survey of dim. Blue Compact Dwarf 
(BCD), and Sm galaxies (Hunter & Elmegreen 2006). They were chosen to be representa- 
tive of the larger dim sample in terms of absolute magnitude My and star formation rate. 
They were also chosen to be relatively free of nearby bright stars that would compromise 
photometry at ultra-low light levels. The new deep imaging data will be referred to here as 
the "deep" data, while the shallower survey data of the larger sample will be referred to as 
"survey" data. The galaxies and several important characteristics are listed in Table [T] 



2. Observations 
2.1. Optical images 

2.1.1. Data acquisition and basic reductions 

The deep optical imaging observations are summarized in Table [2] and K-band and Ha 
images are shown in Figures [l] through Is] The Ha images are part of the survey data (Hunter 
& Elmegreen 2004). Reductions were done using the Image Reduction and Analysis Facility 
(IRAFg. 

In May 2002 we obtained deep B and V images of DDO 133 and V images of IZw 
115 with the Mosaic Imager on the Kitt Peak National Observatory (KPNO) 4-m Mayall 
telescope. The Mosaic Imager consists of 8 2kx4k SITe CCDs and gives a field of view of 36'. 
Because our galaxies are small on the sky, they were centered on one of the CCDs although 
the entire mosaic of CCDs was reduced together. Dome fiats and twilight sky fiats were used 
to determine the pixel-to-pixel variations and scaled to remove the pupil ghost. The images 
were corrected for geometrical distortions, a background offset was removed, and the images 
were stacked, with cosmic ray removal, to produce a single final image in each filter. The 
Mosaic Imager images were calibrated using field images and Landolt (1992) standard stars 
observed with the Lowell Observatory 1.1 m telescope. Because IZw 115 was not observed 
in B, a [B — V)q of 0.4 was assumed to calibrate instrumental V magnitudes. Uncertainties 
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in the calibration are given in Table [2] 

We obtained optical images of DDO 53, DDO 86, and NGC 4163 with the Cassegrain 
Focus Imager on the KPNO 2.1 m telescope. The T2KA CCD was used for V observations 
of DDO 53 and DDO 86, and T2KB was used for observations of DDO 53 in B and of 
NGC 4163 in V and B. The 2048x2048 pixel format yields a 10.4' field of view. Landolt 
(1992) standard stars were also observed for photometric calibration, and the rms of the 
calibrations are given in Table |2} We used observations of an illuminated white spot hanging 
in the dome (dome fiats) to remove pixel-to-pixel variations, and dark sky fiats to remove 
large-scale structure across the CCD. 

Our experience has been that the best method for flat-fielding varies from telescope to 
telescope. Discussion of this point can be found in Massey & Jacoby (1992). In some cases, 
such as at the KPNO 4-m, bright twilight flats can be used to remove both the pixel-to- 
pixel variations and the large scale gradients. But in other cases, such as the KPNO 2.1-m, 
twilight flats do not work well and blank-sky flats work best. Therefore, for the KPNO 2.1-m 
data, we constructed a dark sky flat for each observing run. We obtained multiple images 
with the same exposure time and moon conditions as the galaxy observations. "Blank sky" 
targets were used, but there were still many stars in the images. The telescope was stepped 
between exposures and the series of images was co-added with an algorithm that removed 
the stars by comparing the images pixel by pixel and rejecting pixels that deviated more 
than 3(T from the median, where the a is determined from the noise and gain parameters of 
the CCD. Because the sky flats showed no structure changes from night to night, a single 
sky flat was constructed for an entire observing run. The resulting sky flat was smoothed. 
We normalized the sky flat to a mean value of one near the center, and divided the object 
images by this sky flat. 

In the case of the \^-band image of DDO 53, however, we discovered after the observing 
run that the "blank sky" was contaminated by Milky Way nebular structure (4^ 29™ 45^, 
54°15'33"[J2000]). The Milky Way emission produced diffuse filamentary structure across 
the field of view and so could not be used to construct a dark sky flat. On the other hand, the 
galaxy itself was small and contained in the center of the field of view; even at low light levels 
the galaxy occupied only 14% of the pixels of the image. Therefore, for this galaxy /filter 
only, we used the galaxy images themselves to construct a dark-sky flat: we combined the 
images to produce a single star-free average and fit the non-galaxy portion of the image with 
an order 5 Legendre function, interpolating across the area of the image containing galaxy. 
The galaxy-free region was defined by looking at the extent of the galaxy with the image 
displayed to see faint noise variations in the background. 

We observed the galaxies only at relatively low airmass (< 1.4), high transparency, and 



no moonlight. In all cases we took a series of exposures of the galaxy, and the telescope 
was stepped between exposures in order to move the galaxy around the chip and average 
over flat-fielding imperfections. We only retained images with less than 7% transparency 
differences and with point-spread-function FWHM less than 50% greater than the best in 
the series. The fiat-fielded galaxy images were shifted to align the stars and co-added to 
build signal-to-noise. We used an algorithm to remove cosmic rays and adjusted parameters 
so that the photometric properties, as determined in a simple average of the images, was 
preserved. 

In spite of our great care in determining fiat-field data, the night-sky brightness at 
KPNO (or any dark site) is about 22.5 mag arcsec"^ (see, for example, Neugent & Massey), 
and so achieving 29.5 mag arcsec"^ requires going 7 magnitudes below the sky. No fiat-field 
will be good to 1 part in 6000 (10% photometry) at these light levels. Fortunately, it does not 
have to be. The key to measuring these low light levels is not in the fiat-fielding, but rather 
in the sky subtraction. The fiat-fielding helps in getting the sky relatively flat from one side 
of the frame to the other, but even if there are 1% residual gradients in the flat-fielding, this 
is dealt with by fitting the sky with a 2-dimensional function, as we did. 

Therefore, we edited foreground stars and background galaxies from the images and 
determined a two-dimensional fit of the sky pixels in the edited image — interpolating across 
the galaxy. This procedure works because the galaxies are a small portion of the total field of 
view of the image: approximately 14% for DDO 53, 12% for DDO 86, 10% for DDO 133, 26% 
for NGC 4163, and 4% for IZw 115 and its detached region. It is necessary because there is 
structure in the sky. Typical variations of the background across the image were of order 0.2- 
0.6%. We used a Legendre function with the minimum order that produced a satisfactory 
sky. Typical orders were 3-5. We determined the success of a fit to the sky empirically 
by examining row and column cuts through the sky-subtracted image. We expected the 
background to be flat with noise around zero, and we varied the x and y orders until we 
achieved that. This procedure produced an image of the sky. We then subtracted the sky 
image from the original co-added image and from the edited image to yield sky-subtracted 
images with and without the foreground/background object contamination. 

A legitimate question is whether or not our fits to the sky structure are affected by the 
presence of the galaxy and/or contaminating stars in the field. That can be answered best by 
comparing results obtained on the same galaxy using different telescopes/instruments with 



differing fields of view. This comparison is made in ^2.1.2 We further gain confidence in our 
results in the outer disks of the dwarfs from the fact that we see similar profiles from data 
obtained with different telescopes/instruments and reduced with different types of sky flats 
(DDO 133 and IZw 115 imaged with the Mosaic Camera on the KPNO 4 m and reduced 



with twilight flats, and the other galaxies imaged with the Cassegrain Focus Imager on the 
KPNO 2.1 m telescope and reduced with dark sky flats). 



2.1.2. Surface photometry 

Removal of background galaxies and foreground stars is essential for meaningful azimuthally- 
averaged surface photometry. We attacked this problem in two ways for each image. First, we 
used the edited sky-subtracted images, with some additional masking of small areas in close 
to the galaxy. A problem with this approach is in distinguishing, without radial velocity in- 
formation, stellar clusters in the galaxy from foreground/background object contamination. 
We used the various passband images of the galaxy at our disposal to make the best guess as 
to whether a discrete object in the image was a part of the galaxy or foreground/background 
contamination. The area surrounding the galaxy is a good clue to the type of contamination 
we should expect across the galaxy. Nevertheless, subjective judgement was involved in this 
step. 

The second method we used to account for foreground/background contamination was 
to use the unedited sky-subtracted image and determined an average "sky" value around 
the galaxy that included the foreground/background component. This statistically corrects 
for interlopers, and assumes that the fields around the galaxy represent the contamination 
across the galaxy itself. The problem with this method is that a given annulus in the surface 
photometry might not contain an exactly representative amount of light from extraneous 
objects, and a single foreground star or background object in the faint, outer parts of the 
galaxies can severely affect the surface photometry. For this reason, we prefer the surface 
photometry obtained with the first method of editing contamination from the images, but we 
applied both methods. In fact the two methods produced similar results. The comparison is 
shown in Figures |6] and [7] for the y-band and GALEX NUV photometry, respectively. The 
surface photometry is cut off where there is a large jump in the uncertainty, usually around 
1 magnitude arcsec"^. 

Because dim galaxies are lumpy, we determined the center of the galaxy, position angle 
(P.A.), and ellipticity from an outer isophote in a contour plot of a smoothed version of the 
V image. Then we integrated the photometry in fixed ellipses that increase in semi-major 
axis length. The step size is given in Table [2| and was chosen to be >150 pc and, as seen by 
eye, to cover a sufficient area of the galaxy to average over the most severe local fluctuations. 
We geometrically matched the continuum-subtracted Ha images to the l^-band images and 
applied the same ellipse integrations. The outer most ellipse is shown superposed on the V , 
NUV, and Ha images in Figures [T] through [5] 



We corrected the photometry for reddening using the foreground reddening, given in 
Table III plus 0.05 mag in E(B — V) for internal reddening in B and V and 0.1 mag in 
E{B — V) for internal reddening in Ha. This level of internal reddening is consistent with 
measurements of the Balmer decrement in H ll regions in a sample of 39 dim galaxies (Hunter 
& Hoffman 1999). There, the average reddening in Hll regions is 0.1, and we have taken 
half this to represent the stars outside of Hll regions. We used the Cardelli et al. (1989) 
extinction curve. 

We used a single extinction correction for all radii within each galaxy. This is justified 
by the lack of metallicity gradients in dwarfs (see, for example, Croxall et al. 2009) and the 
correlation between metallicity and attenuation seen in spirals (Boissier et al. 2007). On 
the other hand, little is known about the outer disks of dwarfs, beyond where Ha emission 
is no longer detected, although there are some hints of dust emission even there (see, for 
example, Hinz et al. 2006). Boissier et al., based on /i^ylS" observations, suggest that there is 
no attenuation in the outer disks of two dwarfs (IC1613, WLM), but uncertainties are large 
due to low IRAS fluxes. 

There is a region of stellar emission near, but detached from, IZw 115 that we did not 
include in the galactic photometry. It is located at 15^32™51.0" 46°25'51" (2000), 101.5" 
(=6.0 kpc) from the center of IZw 115, southwest along the galaxy's major axis, and is 
marked in images in Figure [S] In a circular aperture of radius 14.85", the region has a V 
magnitude of 20.39±0.04 and a surface brightness fiy of 26.9 mag of one arcsec^. If it is at 
the distance of IZw 115, it has an My of —10 and a diameter of several kpc. Thus, it could 
be a very small dwarf companion. 

DDO 133 is clearly barred, and the outer isophote that was used to determine the 
galaxy morphology is located beyond the bar. The bar has a P.A. that is —13°, an 11° twist 
from the outer isophote. The center of the bar is located at 12^32'"54.5" 31°32'27", 18.8" 
(=560 pc) from the center of the outer isophote. The minor-to-major axis ratio b/a of the 
bar is 0.44, compared to 0.69 for the outer isophote, so the outer isophote is rounder. The 
major axis of the bar is 133" (=3.9 kpc) long. We experimented with changing the centering 
and P.A. of the photometry ellipses, but found that there was no significant change in the 
photometry for a reasonable alternate choice of the center (10" offset to the southwest) and 
P.A. (5° counter-clockwise). 

DDO 86 presents a special challenge to determining the position angle of the major axis. 
At an outer isophote [nv ~ 27.2 mag arcsec"^), the P.A. of the major axis is 85° and the 
minor-to-major axis ratio b/a is 0.86. This is close to the galaxy morphology determined 
in the survey data: P.A.~72° and b/a ~0.82 (Hunter & Elmegreen 2006). However, in the 
deeper imaging data that we have here, we see that there is a faint extension of starlight 



to the north. At an isophote level that includes this extension (/ly ~ 28.4 mag arcsec"^), 
the P.A. of the galaxy rotates clockwise to 32.5°, although the shape remains the same 
{b/a ~ 0.87). An inner isophote (/iv ~ 25.6 mag arcsec"^) has a major axis that is rotated 
counter-clockwise from both of these (113.5°) and the galaxy is more elongated {h/a ~ 0.65). 
Finally, the kinematical major axis determined from the Hi velocity field has a P.A. of 
—46°. It is possible that this sequence of optical P.A.s — rotating clockwise from higher 
surface brightness to fainter — and change in shape — becoming rounder in the middle and 
outer parts — represents an offset bar; the center of the inner isophote is located 360 pc to 
the southwest of the center of the extreme outer isophote. For the surface photometry we 
have chosen to use the P.A. from our H I kinematic analysis. 

Since we have surface photometry from the survey images of each of these galaxies, 
we compare the deep imaging photometry with that from the survey data in Figure [6j We 
see that the agreement is quite good in most cases. In addition, the breaks in the surface 
brightness profiles seen in the survey data of DDO 86 and DDO 133 are there at the same 
radius and surface brightness level in the deep data. And, the azimuthally-averaged B — V 
colors in DDO 133 match closely the survey data in the region of overlap. The exception to 
the good agreement is DDO 53. Beginning at a radius of 1.2' (1.25 kpc) and a /xy of 25.8 mag 
arcsec"^ the new surface photometry is systematically fainter. The differences are greater 
than the uncertainties in the photometry. The survey B — V color is also systematically 
redder than the deep photometry beginning at a radius of 0.77' (0.81 kpc; /iy ~ 24.9 mag 
arcsec"^). In the survey data, there is a gradient in B — V , with the galaxy becoming redder 
with radius. In the deep data, the color is relatively constant. We note that the integrated 
B — V oi DDO 53 from the deep image (0.30±0.01, not corrected for reddening) is close 
to that given by de Vaucouleurs et al. (1991; 0.35), while that from the survey data (0.48) 
is significantly redder. Therefore, we suspect that the deep photometry is better than that 
measured from the survey data for this galaxy in the region of overlap. 

From the ellipse photometry we can derive several quantities. First, the last ellipse yields 
integrated photometry and that is given in Table |3} Second, fits to the surface brightness 
profiles yield disk scale lengths Rd and central surface brightnesses /xy^o- These are listed 
in Table |4| Four of the 5 galaxies— DDO 53, DDO 86, DDO 133, and NGC 4163— exhibit 
double exponentials. The profile in the outer disk falls off more steeply than that in the 
inner disk in DDO 53, DDO 86, and DDO 133. In NGC 4163 the outer profile is shallower 
than the inner. DDO 86 and DDO 133 were already known to have double exponentials from 
the survey data, but the survey data of DDO 53 did not go out far enough to give a hint 
of the break in the profile. The break in NGC 4163 is more subtle than in the other three 
galaxies. Nevertheless, the scale lengths are different — 0.26 ± 0.02 kpc versus 0.31 ± 0.004 
kpc, and the central surface brightnesses implied by each differ by 0.54 magnitude. For those 
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four galaxies with breaks, the Rd and ixyfi are given for both the inner and outer parts of 
the profiles, along with the break radius Rbt at which the slope changes. In addition DDO 
53 and IZw 115 show breaks at small radii, a little over one disk-scale length that coincide 
with concentrations of star formation. These are marked in the surface brightness plots, but 
not included in Table |4j The surface photometry and color profiles are shown in Figures |8] 
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2.2. Ultraviolet images 

We obtained deep images of our sample galaxies with the Galaxy Evolution Explorer 
satellite {GALEX; Martin et al. 2005) to trace and characterize star formation into the outer 
disks of dwarf galaxies where Ha may not be an effective tracer of recent star formation 
(see, for example, Thilker et al. 2005, 2007a; Boissier et al. 2007). The UV — ultraviolet 
light coming directly from OB stars — can easily trace young stars in the outer galaxy, and 
the UV has the advantage in a patchy star-forming environment of integrating over a longer 
timescale than does Ha (10 Myrs). 

The galaxies are listed in Table |5] along with the exposure times and tile name of the 
GALEX images. GALEX imaged simultaneously in two channels: FUV with a bandpass 
of 1350-1750 A, an effective wavelength of 1516 A, and a resolution of 4.0" and NUV with 
a bandpass of 1750-2800 A, an effective wavelength of 2267 A, and a resolution of 5.6". 
The images were processed through the GALEX GR4/5 pipeline and were retrieved as final 
intensity maps with a 1.5" pixel scale. The GALEX Held of view is a circle with 1.2° diameter, 
and we have extracted a portion around our target galaxies. 

We corrected the UV photometry for extinction using the same total reddening E{B—V)t 
as was used for correcting the optical. We combined the E{B — V)t with the extinction law 
of Cardelli et al. (1989) to produce the extinction Apuy = 8.24E(i? — l^)^, and, interpolating 
over the 2175 A bump, A^uv = 7.39E(i? — V)t- However, the FUV extinction law is known 
to vary from galaxy to galaxy, and the NUV filter straddles the 2175 A extinction feature, 
which also varies not only from galaxy to galaxy but from place to place within galaxies 
(for example, Gordon et al. 2003). Star formation activity appears to play as important a 
role, perhaps even more important, than metallicity in determining the UV extinction curve 
and the strength of the 2175 A bump. In the metal-poor SMC, for example, Gordon & 
Clayton (1998) have shown that lines of sight through the actively star forming bar have 
UV extinction curves that are linear with 1/A, rising more steeply into the FUV than in 
the Milky Way, and no 2175 A bump is present, while a sight line through the quiescent 
wing has an extinction curve that is similar to that in the Milky Way with a shallower 
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slope into the FUV and a 2175 A bump. In our dwarf galaxies, the level of star formation 
activity is highly variable from galaxy to galaxy and from place to place within a galaxy, 
potentially resulting in a variable FUV extinction law. Wyder et al. (2007) have dealt with 
this for a large sample of galaxies by adopting the Cardelli et al. extinction law for the 
Milky Way and determining the best average reddening from convolving theoretical galaxy 
spectral energy distributions (SEDs) with the GALEX Glter transmission curves. In Hunter, 
Elmegreen, & Ludka (2010), we examined the results of convolving different SMC, LMC, 
and Milky Way extinction curves (Cardelli et al.l989, Gordon et al. 2003) with constant 
star formation rate galaxy SEDs constructed from the Bruzual & Chariot (2003) stellar 
population library and with the GALEX filter transmission curves. The Ax/'E{B — V) 
varies by 1.5 magnitudes depending on the extinction curve. Considering the variety of 
environments represented by our sample and the fact that extinction is minimal except for a 
few galaxies with higher than average foreground extinction, we decided to adopt the Wyder 
et al. extinctions: Apuv = 8.24E(_B — V)t and Aj^uv = 8.2E(_B — V)t. For most of our 
galaxies, the variations represented by alternate extinction laws result in differences of order 
0.05 mag. 

Surface photometry of the NUV and FUV images proceeded in the same manner as 
for the optical images. We edited foreground/background objects from the images, fit the 
cleared area around the galaxy, and subtracted this from the edited and from the unedited 
images. For the unedited sky-subtracted image, we determined the average background 
light from contaminating objects and subtracted this from the image. We then geometrically 
transformed the UV images to match the V^-band scale and orientation and measured surface 
photometry on the edited and unedited images using the same ellipse parameters as for the 
V^-band images. The two methods of accounting for background/foreground contamination 
are shown in Figure [7j The results of the two methods agree well. Integrated photometry are 
given in Table |3j The surface photometry and color profiles are shown in Figures |8] through 

m 

Since we are looking at faint light levels in outer disks, we need to consider whether the 
UV light we see truly comes from young stars or whether there might be other sources of 
light that dominate out there. One such potential source is white dwarf stars. DDO 133, 
as a test case, has an FUV — NUV color of order 0.1 mag in the outer disk, not counting 
the colors with relatively high uncertainties. A white dwarf of that color has an effective 
temperature of order 16,000 K. At a distance of 10 pc, a star with logg = 6 and a radius 
R = Rq would have an absolute NUV AB magnitude of Mj^my = 1.84 (L. Bianchi, private 
communication). Assuming a mass for the white dwarf of 0.7 Mq , the star should have a 
radius of O.IRq . Thus, the apparent NUV magnitude would be 35.7. To reproduce a fiNuv 
of 28 mag arcsec"^ at a radius of 2i?£) in DDO 133, we would need about 1.4 white dwarfs 
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pc~^. That is 2 times higher than the density of white dwarfs in the solar neighborhood. 
For white dwarfs with smaller radii or higher logg, the number of white dwarfs that are 
needed goes up, so this is probably a lower limit. However, the solar neighborhood should 
have a much higher stellar density than the outer disk of a dwarf. In the solar neighborhood 
Hv is 22.7±0.2 (Ishida & Mikami 1982), and that is a factor of 30 higher than the surface 
brightness in DDO 133 at the radius we are considering (/zy = 26.5). So, the stellar density 
of white dwarfs needed to account for the UV light in DDO 133 fails by a factor of about 60. 
Even 3 magnitudes fainter — at the edge of our detection, the numbers needed are 4 times 
too high and so it is unlikely that white dwarfs alone could produce the UV starlight that 
we see. We also considered horizontal branch stars since these dominate the UV light in, for 
example, M32 (Brown et al. 2000). However, the FUV — NUV color of these stars are of 
order 2 and redder. This color is far redder than what we see in dwarfs. 



2.3. Hi 

Three of the galaxies in our sample — DDO 53, DDO 133, NGC 4163 — are part of the 
LITTLE THINGS {LITTLE: Local Irregulars That Trace Luminosity Extremes; THINGS: 
The HI Nearby Galaxy Survey) project, a multi-wavelength survey of dwarf galaxies (LITTLE 
THINGS: Hunter et al., in preparation; THINGS: Walter et al. 2008). LITTLE THINGS 
has H I observations at additional VLA arrays, including B array, that provide higher spatial 
resolution. They will provide much better maps for analyzing the relationship between the 
gas and star formation. Here we restrict ourselves to the data obtained for this project and 
examine general correlations between the stars and the gas. 



2.3.1. Calibration and Mapping 



We obtained Very Large Array (VLA** ) observations of H I -line emission in our galaxies. 
The dates and configurations of the L-band observations are given in Table [6j The total 
bandwidth was 1.56 MHz, and 2 IFs were used. The data were on-line Hanning smoothed, 
and the channel separation was 12.2 kHz (2.6 km s^^). The data from 2006 and 2008 were 
obtained when the VLA was in transition to the EVLA. Unfortunately, in converting the 
signal from the EVLA antennas, power was aliased into the lower 500 kHz of the bandwidth. 



** The VLA is a facility of the National Radio Astronomy Observatory (NRAO). The National Radio As- 
tronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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The consequence of this is that all EVLA-EVLA baselines had to be flagged and not used 
(see Ficut-Vicas et al., in preparation, for a discussion of how to deal with this problem). 

We subtracted the continuum emission in the uv-p\ane using line-free channels on either 
end of the spectrum and, in the case of DDO 53, combined the data from the two arrays. 
To make maps, we employed a routine in NRAO's Astronomical Image Processing System 
(AIPS) that enables one to choose a sample weighting that is in between the standard 
"natural" weighting, which gives the highest signal-to-noise but at the expense of long wings 
in the beam, and "uniform" weighting, which gives the best resolution but at the cost of 
signal-to-noise and the presence of negative sidelobes. We chose a sample weighting that gives 
a formal increase of only 5% in noise yet with a significant improvement in beam profile over 
what would have been achieved with "natural" weighting. The resulting synthesized beam 
FWHM are given in Table |6| We deconvolved the maps until there were roughly comparable 
numbers of positive and negative components. 

To remove the portions of each channel map without emission that are only contributing 
noise to the map, we used the maps smoothed to twice the beamsize for a conditional transfer 
of data in the unsmoothed maps. The channel maps were blanked wherever the flux in the 
smoothed maps fell below 2.5a. Flux-weighted moment maps were made from the resulting 
data cube. Integrated emission — moment zero — maps are shown in Figures [13] through 



17, along with velocity field (moment one) and velocity dispersion (moment two) maps. 
Integrated H I contours superposed on the \^-band images are shown in Figures [T] through 

El 



2.3.2. Velocity fields 
Position velocity diagrams (panel e), velocity fields (panels b and d), and velocity disper- 



sions (panel c) are shown in Figures 13 to 17 Panel (b) shows the standard intensity- weighted 
mean velocity field, and panel (d) shows the velocity field that results from fitting profiles 
with a Gauss-Hermite hs polynomial. 

Because of the possible presence of multiple velocity components, we have deconvolved 
ordered from non-ordered motions using the method described in Oh et al. (2008). First, 
an approximate rotation curve is derived from the major axis position-velocity diagram and 
this is used to create an artificial velocity field. Then in an automated method, we fit single 
Gaussians to the major velocity component in each profile across the galaxy. The velocities 
of these components are compared to the artificial velocity field. If they deviate more than 
a prescribed amount, the program examines the possibility that the bulk motion at that 
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point is represented by a secondary velocity component. From this, we extract the bulk 
velocity field that then becomes the input velocity field for the next iteration, which is now 
fit using the "tilted ring model" rather than the major axis position-velocity diagram. This 
loop continues until the mean difference between successive velocity fields is less than three 
times the channel width. This gives a final bulk velocity field. Parameters derived from this 
process (center position, P. A., inclinations) are summarized in Table [TJ Following Bureau 
& Carignan (2002) we corrected the rotation velocities for asymmetric drift, and the final 



rotation curves are shown in the graphs in Figures p^ through 12 



The velocity fields of DDO 86 and DDO 133 exhibit well-behaved ordered rotation. 
In IZw 115 the velocity field is less well-behaved, but still ordered motion is present. The 
velocity fields in DDO 53 and NC 4163, on the other hand, are quite peculiar. Here we 
discuss each galaxy individually: 

DDO 53. DDO 53 is a part of the THINGS study of H I in nearby galaxies (Walter 
2008) and has been analyzed by Oh et al. (in preparation). The kinematical axis of DDO 53 
is at a position angle of 155°, while the optical major axis is at an angle of 81°, a difference 
of 74°. Furthermore, the fit to the Hi velocity field yields a variable inclination angle. The 
kinematic center and the optical center are 32" apart, with the kinematic center being NE of 
the optical center. The fit to the velocity field is quite poor, and the rotation curve is only 
meaningful out to a radius of 0.9 kpc. 

DDO 86. The bulk velocity field of DDO 86 is well determined with clean fits to the 
position angle of the kinematic axis, the center position, and the inclination of the H I disk. 



The optical, however, is more problematic, as discussed in ^ 2.1.2 The V^-band major axis 
rotates clockwise from higher surface brightness to fainter. This may indicate the presence 
of an off-center bar. So, for the optical and UV surface photometry we used the P.A. derived 
from the H I kinematics. The optical center determined from outer isophotes agrees within 
14" with the H I center, with the H I center being NW of the optical. 

DDO 133. The bulk velocity field of DDO 133 is well-behaved. The kinematic P.A., 
inclination, and center agree very well with those determined for the optical. 

NGC 4163. Although rotation is seen in the gas, the velocity field is so peculiar that 
parameters for the tilted ring model — center, P.A., and inclination-were taken from the 
optical. 

IZw 115. IZw 115 shows rotation, although the velocity field is rather peculiar. The 
position angle of the kinematic axis was found to vary between 50° and 78°, while the optical 
major axis is at an angle of 41°. On the other hand, the inclinations agree within 7°, and 
the centers, within 4". The detached region seen in the optical and NUV (Figure [s]) shows a 
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hint of rotating independently of the main body of IZw 115. Panel (d) in Figure 17, which 
shows the velocity field resulting from fitting profiles with a Gauss-Hermite h^ polynomial, 
can be interpreted as showing overall rotation with an axis at a position angle of roughly 
90°, clearly distinct from IZw 115. This could indicate that the detached region is a small 
companion galaxy that is interacting with IZw 115. 



2.3.3. Surface density profiles 

The azimuthally-averaged H I surface densities T^m and velocity dispersions a-^ei of the 
rotational component were determined using the center, P. A., and inclination angle given in 
Table [7j The T^hi was corrected for He to produce the gas surface density Tigas- The gas 
surface density profiles and average velocity dispersions are shown as a function of radius in 



Figures pi through 12 



3. Stars in the outer disks 

3.1. Where do dim stellar disks end? 

We have traced starlight in the \^-band to ultra-low surface brightness levels. Yet, 
the stellar surface brightnesses in the outer disks of our 5 galaxies continue to decline ex- 
ponentially as far as we can measure it. That is, we run out of signal-to-noise before we 
convincingly run out of galactic starlight. We have not found the edges of the stellar disks. 
This is consistent with results from stellar photometry that have traced the stellar disks in 
the LMC (Saha et al. 2010), M33 (Barker et al. 2011), and M31 (Richardson et al. 2008) to 
even more extreme radii. 



3.2. Breaks in the Surface Brightness Profiles 

The surface photometry of DDO 133 shows a change in the slope at 1.8-R/). Beyond 
there, /xy drops more steeply. Although less extreme, there are noticeable changes in the 
slopes in the outer disks of DDO 53, DDO 86, and NGC 4163 as well. And, in DDO 86 and 
DDO 133 the breaks are also apparent in NUV surface brightness profiles. The break in NGC 
4163 is coincident with the outer limit to a central bump in NUV and Ha emission. And 
in DDO 86 and DDO 133, the outer disk breaks occur where Ha ends or declines sharply. 
In addition to the breaks at several disk scale lengths, there are breaks or variations in the 
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surface brightness profiles at smaller radii. We see these in DDO 53 and IZw 115. The 
inner bump in the l^-band surface photometry in DDO 53 is mimicked in NUV and Ha and 
corresponds, more or less, to a change in FUV— NUV and B — V colors. The flattening out 
of the \^-band surface brightness profile in the central 1.5Rd in IZw 115 is accompanied by 
a flattening in the NUV proflle, and the location of a single bright Hll region. The only 
break that does not coincide with a change in Ha is the outer break in DDO 53. Thus, most 
of these breaks appear to be connected with changes in the star formation activity. 

Double-exponential surface brightness proflles are also seen in about 1/4 of our survey 
of dim and Sm galaxies (Hunter & Elmegreen 2006) and in the outer parts of many spirals 
(van der Kruit & Shostak 1982; Shostak & Van der Kruit 1984; de Grijs et al. 2001; Kregel 
et al. 2002; Pohlen et al. 2002; Kregel & van der Kruit 2004), including one low luminosity 
spiral (Simon et al. 2003). It is now also seen in disks at high redshifts (0.6 < z < 1.0; 
Perez 2004). Traditionally, these surface brightness breaks have been called "truncations," 
but this implies that the break signals the edge of the stellar disk. In dwarf galaxies, this 
does not seem to be the case; the break is a rather abrupt change rather than an end. But 
what changes at that radius? Elmegreen & Hunter (2006) have suggested that the primary 
driver of star formation changes at the break, from gravitational instabilities in the interior to 
turbulence compression and supernovae in the outer disk. The nature of the breaks in dwarfs 
is explored further with our much larger survey data by Herrmann et al. (in preparation). 



3.3. What is the star formation history in the outer disks? 

What does our observational limit of 29.5 mag arcsec"^ in fiy mean? Let us take DDO 
133 at a distance of 6.1 Mpc as an example. An annulus centered at /xy = 29.5 mag arcsec"^ 
with a width of 1 kpc would have an area of 45617 arcsec^. This gives an My of —11.1 
and an Ly of 2.3 x W^Lq . We measure a B - V of order 0.3-0.4 in DDO 133. A B -V 
of 0.35 gives a mass to light ratio M/Ly of 1.07, under the assumption of a constant star 
formation rate and Salpeter (1955) stellar initial mass function (Bell & de Jong 2001). Using 
this stellar mass to light ratio, we obtain a mass of stars in the annulus of 2.5 x 10® M0 . 
The stellar density is 0.06 M0 pc~^, or one Sun every 16 pc^. 

To look at ages of the stars in the outer disks, we consider an FUV— NUV color of 0.1, 
which is the value measured in the outer disk of DDO 133 (1-1.7/?/)), before the uncertainties 
in the color rise. At solar metallicity, an FUV— NUV color of 0.1 magnitude, implies a single 
age population that is 300 Myrs old (Bruzual & Chariot 2003). But at Z = 0.004 the age 
increases to 670 Myrs and at Z = 0.0004 to 1 Gyr. At lower metallicities, stars become bluer 
in FUV— NUV and so older stars take on the colors of that one would flnd for younger stars 



17 



in a higher metalhcity population. Thus, at our metalhcities, the stellar populations are not 
likely to be as young as we might otherwise predict. 

To look at the rates necessary to form the stars we find in the outer disks, we have 
computed azimuthally-averaged profiles of star formation rates (SFRs) derived from Ha, 
FUV, and V-band luminosities: SFRua, SFRpuy, and SFRy. We use the formulae outhned 
in Hunter et al. (2010) with a Salpeter (1955) stellar initial mass function (IMF) and assuming 
a constant SFR over long timescales. For the FUV and Ha formulae, we use the formulae 
from Kennicutt (1998) modified for the sub-solar {Z = 0.008) metalhcities appropriate for 
dwarf galaxies. To determine the SFR averaged over a longer time scale, we use the V- 
band luminosity to determine the mass in stars and assume a timescale of 12 Gyr for the 
formation of that stellar mass. We use a stellar M/Ly ratio from Bell & de Jong (2001) that 
is a function of the {B — V)o color. If there is no {B — V)o measurement or the uncertainty 
in the color is > 0.2 magnitude, we assumed the last good measurement oi B — V. In the 
case of DDO 86 and IZw 115, for which we do not have deep 5-band images, we used the 
survey data. For DDO 86 -B — V is constant with radius in the survey photometry; for IZw 
115 the color reddens with radius and we used the survey photometry as far out as it goes. 
The SFR determined from the l^-band is the SFR necessary to produce the original mass 
in stars. We correct this by a factor of 2 for the mass recycled back into the interstellar 
medium from Brinchmann et al. (2004). 

Now we consider the effects of metalhcity and a further correction to the SFR formulae 
based on the actual metalhcities of our sample. Two of the galaxies in our sample have 
oxygen abundances 12 + log 0/H measured from H ll regions: The measured value of the 
oxygen abundance in DDO 53 is 7.5±0.1 according to Hidalgo-Gamez et al. (2003) and 
7.82 ± 0.09 according to Croxall et al. (2009). Given the coursness of the metalhcity grid of 
the models, both measurements give the same SFRs. In DDO 133 the oxygen abundance is 
8.0±0.05 (Hunter & Hoffman 1999). For the other galaxies, there are no measured values, so 
we have estimated the oxygen abundance from Mb and Richer & McCall's (1995) correlation 
between Mb and 12 + log 0/H. This yields 12 + log 0/H of 8.0 for DDO 86 and IZw 115 and 
7.7 for NGC 4163. Metalhcities are unknown in the outer disks, but are likely to be lower 
than these values. To modify our star formation formulae for these metalhcities, we use the 
synthesis models of STARBURST99 (Leitherer et al. 1999), with constant star formation 
rates and Salpeter stellar IMFs. Their Figure 54 gives the luminosity at 1500 A for Z of 
0.008, 0.004, and 0.001, and this allows us to adjust our Z = 0.008 formulae for FUV-based 
SFRs to lower metalhcities. Similarly, their Figure 78, which gives the number of photons 
below 912 A, enables us to correct our Ha formula, and finally Figure 58 allows us to correct 
B — V, which enters into the formula for M/Ly and the \/-band SFR. 



For our galaxies, the SFRs at a radius where nv = 29.5 mag arcsec"^ are given in Table 
|8] Only SFRpuY and SFRy are given since Ha is not detected at those radii. In DDO 133, 
for example, the stars continue for at least another 1.5 kpc beyond where Ha emission is 
no longer detectable. In DDO 133 at an annulus centered at ^v = 29.5 mag arcsec"^, the 
SFR determined from the stellar mass SFRy is 0.00001 Mq yr~^ kpc~^. That corresponds 
to 0.0004 Mq yr~^ in the 1-kpc wide annulus we are using as an example. This corresponds 
to approximately 6 Orion Nebulae over an area of 40 kpc^. This is the level of activity that 
we should expect to find on average at any given moment. 



The azimuthally-averaged SFR profiles are shown in Figure [18] as ratios of SFR}ia and 
SFRpuY to SFRy. We are able to measure SFRs from FUV and l^-band images out to 
3-7 disk scale lengths. Again, Ha is not detected as far out as V and FUV emission, and 
so is not a good indicator of star formation in the outer disk (see Hunter et al. 2010 for a 
discussion of this issue, as well as Meurer et al. 2009 and Lee et al. 2009). However, in the 
region where Ha and FUV emission overlap, the two measures of SFRs track each other well 



except in NGC 4163. This is obvious to the eye as parallel curves in Figure 18 In addition, 
within the likely systematic uncertainties, the absolute levels of the Ha and FUV SFRs are 
the same except in the center of DDO 53, within 1.5Rd of NGC 4163, and in IZw 115. In 
those cases the Ha SFR is significantly lower than the FUV SFR, implying a higher SFR 
several 100 Myrs ago compared to the most recent 10 Myrs. 

These three measures of the SFR can track star formation over different time scales: 10 
Myrs, several times 100 Myrs, and over 1000 Myrs, respectively for SFRjia, SFRp^y, and 
SFRy. What then do these profiles tell us about the radial variations in the star formation 
histories in our five galaxies? 

DDO 53. In the central 2Rd of DDO 53 the SFR over the recent 10-300 Myrs has 
been higher than the average in the past by factors of 3-4, while in the outer disk the SFR is 
suppressed by similar factors compared to the past average. Recent star formation in DDO 53 
has retreated to the central regions. This is also reflected in the color profiles: FUV— NUV is 
constant in the inner disk scale length and becomes progressively redder beyond that; B — V 
is slightly bluer in the central regions. 

DDO 86. Star formation in DDO 86 has been approximately constant over the past 
Gyr at all radii with only a slight depression of the current SFR relative to the average past 
rate in the center. Similarly, the FUV— NUV color is nearly constant with radius. 

DDO 133. In the central Q.hRo of DDO 133 the SFR over the recent 10-300 Myrs is 
depressed compared to the average in the past by factors of several, but beyond that the 
SFR has been roughly constant until it drops a bit beyond 2i?£). The FUV— NUV color is 
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red in the center and becomes steadily bluer with radius until it turns a bit redder again at 
2Rd. B — V, on the other hand, is nearly constant with radius to 2.5Rr,. This flat color 
profile implies that the longer timescale star formation history has been the same throughout 
the entire galaxy, including the outer edges. 

NGC 4163. In the central IRd the SFR averaged over the past several 100 Myrs is 
comparable to that measured over the past 1 Gyr, but the SFR measured over the past 10 
Myrs is down by a factor of about 10. Thus, in the center of NGC 4163 star formation 
has recently declined significantly. Outside of the central IRd, the recent star formation 
activity is depressed relative to the past average by a factor of up to 30. It seems that only 
the center of this galaxy is currently engaged in significant star formation. In terms of the 
colors, FUV— NUV becomes steadily redder with radius beyond IRd- B — V becomes redder 
to 2.5Rd and then is constant with radius. 

IZw 115. IZw 115 is a Blue Compact Dwarf, which are often assumed to be starburst 
or post- starburst systems. In the case of IZw 115, however, the SFRpuv is comparable to 
the average past rate in the central 2Rd and is depressed by factors of a few in the outer 
disk. However, the most recent star formation, as seen in Ha, is down by a factor of order 10 
where it is detected in the inner 2.5Rd. Thus, star formation in the center of IZw 115 was 
much higher in the past several 100 Myrs than it is now. In that sense, it is post-starburst. 
The FUV— NUV color refiects this, becoming steadily redder with radius. 



4. Gas Morphology and Kinematics: No Normal Galaxies 

4.1. DDO 86 and a Hole in the ISM 

The Hi distribution in DDO 86 is dominated by a central hole with a higher density 
gas ring around it. The hole is not as distinct as a simple circular region of low gas density; 
there are several minima and some gas "clouds" within the shell. But the general hole/ring 
morphology is clearly present. The hole, measured from the inner edge of the shell, has a 
diameter of order 5.2 kpc, the structure, measured from the midpoint of the shell, has a 
diameter of order 6.9 kpc, and the shell itself is of order 2.0 kpc thick. The mass of Hi 
in the shell, about 4 x 10^ Mq , is 41% of the total mass of the galaxy. Many other dwarf 
galaxies exhibit this type of morphology (see, for example, Dopita et al. 1985; Puche et al. 
1992; Martin 1998; Walter & Brinks 1999; Parodi & Binggeh 2003; Cannon et al. 2005; de 
Blok & Walter 2006; Cannon et al. 2011a,b). One galaxy with a similar appearance and 
similar hole/shell dimensions and relative mass is DDO 88 (Simpson et al. 2005). 

Although DDO 86 has a beautifully normal bulk rotation velocity field, at almost every 
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point within the shell and hole there are weak non-circular motions as well. In particular 
there appears consistently to be a component that is ~10-25 km s~^ off the rotation speed at 
a given position. This suggests that the shell is young enough to still be expanding at 10-25 
km s~^. McCray and Kafatos' (1987) model of stellar blown shells give an age of 60-150 
Myrs for a shell observed with an expansion velocity of 10-25 km s~^ and a radius of 2.6 kpc. 
If the ambient gas density before the shell formed was of order 0.1 cm^^, 100-1600 stars with 
masses greater than 18 Mq would have been required to produce the shell. That means 
that a substantive burst of star formation took place in the center of DDO 86 in the past 
200 Myrs. 



IS con- 



However, neither the FUV— NUV color in Figure [9] nor the SFRs in Figure 18 
sistent with this picture. The Hi minima regions in the center yield an FUV— NUV of 
0.24±0.02, which indicates an age of order 550 Myrs for a Chabrier stellar initial mass func- 
tion, a single starburst population, and a metallicity of 0.008 (Bruzual & Chariot 2003). 
A single age stellar population at 200 Myrs of age should have an FUV— NUV color of or- 



der —0.4. Similarly, the SFR profiles shown in Figure 18 are counter to expectations: The 
SFRiT'f/v/SFRy and SFRji/^/SFRy ratios should show a bump up in the center, indicating a 
higher recent SFR compared to the average over a long timescale. Instead these ratios show 
a slight dip in the center. These discrepancies could be due to a combination of factors: 
for example, a more complex star formation history (see, for example. Cannon et al. 2011), 
contamination of the younger stellar population by an older underlying disk, local extinction, 
and UV calibration uncertainties (see Zhang et al., in preparation). 



4.2. Stellar Bars 

One of our galaxies clearly has a stellar bar and a second galaxy probably has one: The 
bar in DDO 133 is obvious from a casual look at the \^-band image, while that in DDO 86 
is suspected because the deep \^-band image shows a different isophote shape compared to 
shallower images. Both of these galaxies have normal bulk rotation fields. In DDO 133 we 



see a slight S-shape to the iso- velocity contours (Figure 15) that is characteristic of a barred 
disk. Such a feature is not obvious for DDO 86 (Figure 14). Interestingly, DDO 133 also has 
UV-bright clumps all the way out to 2 disk scale-lengths and large UV-bright regions are 
found especially at the ends of the bar. H I profiles are double-peaked in the H I minimum 
that corresponds to a concentration of UV clumps on the northeastern edge of the bar. And 
Hi profiles show broad features, probably due to multiple components, in the Hi minima 
to the east of the bar and in the region of bright UV emission to the south of the bar. 
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4.3. Peculiar Kinematics 

Two of our galaxies are particularly kinematically peculiar. Extracting the bulk rota- 
tion from DDO 53 and NGC 4163 was difficult and the rotation curve for DDO 53 is not 
meaningful outside a radius of 0.9 kpc (but see Oh et al. 2011). Furthermore, the rotation 
speeds that were measured are small, less than 10 km s~^. In the outer disk, the Hi as- 
sociated with NGC 4163 sports two "spiral arm" like features — one on the southern edge 
extending to the west and the other on the northern edge pointing to the NE (Figure |4]). 



In addition, the iso- velocity contours show a strange extension to the SW (Figure 16). Hi 
profiles show complex, multiple (3-4) components outside of the central H I bright region, 
especially in the H I extension to the South of center and along the ridge to the West of the 
central region. In DDO 53, the region of highest velocity dispersion is found along a grand 
arc on the SE edge of the H I distribution, where the iso- velocity contours are most bunched 



together (Figure 13). Except for the southern of the two Hi clumps, the Hi profiles are 
everywhere broad and complex with multiple components. In neither galaxy is it obvious 
from the optical morphology that anything unusual is going on in these galaxies other than 
the fact that the current star formation is concentrated to the center in both systems. 

IZw 115 potentially has an excuse at disturbed kinematics, and the bulk velocity field, 
although more obvious than in DDO 53 and NGC 4163, is not clean. The detached region to 
the SW (Figure [5]) could be a tiny companion interacting with IZw 115. The gas associated 
with IZw 115 extends to that region, and the region itself seems to be rotating at a different 
position angle from the main body of IZw 115. The signal-to-noise is low, but the Hi 
profiles are everywhere complex, indicating multiple (3-4) velocity components. IZw 115 has 
one central star forming complex, as is characteristic of Blue Compact Dwarfs. 



5. Comparison of stars, gas, and star formation 



To compare the stars and gas, in Figures [8] through [12] we show the azimuthally-averaged 
surface photometry plotted against radius from the galaxy center. The top panel in the graph 
includes properties of the gas: the gas surface density E^^s and velocity dispersion of the 
rotational component a^ei- The second panel shows the optical and ultraviolet properties: 
V, H«, and NUV surface photometry. The third and fourth panels show optical and UV 
colors: {B — V)o and {FUV — NUV)o- From comparison of Hi and optical images (Figures 
[1] through Figure Is]), we see that the Hi does not extend much further than the stars in 
DDO 53 and DDO 133. In DDO 133, there are bright UV regions that extend almost to 
the edge of the observed H I , and the big UV-bright complexes are found roughly near the 
ends of the stellar bar. In DDO 53, NGC 4163, and IZw 115 the recent star formation is 
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concentrated to the center, and we do not find any UV-bright knots in the outer disk. In 
DDO 86 there is a central hole in the gas. Many of the UV knots are associated with clouds 
around this hole. 



In Figure 19 we plot the azimuthally-averaged \ogT,gas against logSFR for the three 
star formation measures: SFRjT'fjy, SFRj/q, and SFRy. This is similar to plots by Bigiel et 
al. (2010) for pixel-by-pixel distributions. We see a general trend of higher SFR for higher 
gas density except for some radii in DDO 53 and NGC 4163. Those two galaxies are also 
are most peculiar kinematically. 

The Toomre (1964) model of a differentially rotating thin disk appears to explain star 
formation in the inner disks of spirals (Kennicutt 1989). In this model there is a critical gas 
density Sc above which the disk is unstable to ring-like perturbations in the radial direction. 
Gas that is above Ec will form star-forming clouds; gas below Ec is stable to cloud formation. 
The ratio of S^ to Sc is shown as a function of radius in the five galaxies in the top panels 



of the graphs in Figures p^ through 12 S^/Sc is everywhere less than one. The ratio is 
extremely low throughout DDO 53 (<0.03), but that may be due to the extremely low 
rotation speed in that galaxy. The ratio of S^/Sc is low (<0.3), but roughly constant where 
it was measured over the disk of IZw 115. For the other three galaxies — DDO 86, DDO 133, 
and NGC 4163 — S^/Sc rises to a peak outside the center and then declines to the outer disk. 
In DDO 133, for example, the NUV surface photometry extends to the regime where S^/Sc 
is 0.05, well below the value to permit star formation. So, how is star formation proceeding 
at very low gas column densities? 

One possibility is that clouds of high enough density for star formation form in outer 
gas disks through turbulence compression (Elmegreen & Hunter 2006). Our Hi maps are of 
low spatial resolution, but along the northern edge of DDO 86, for example, there appear 
to be two Hi clouds that stand out. Both have peaks that are 6a above the surrounding 
gas. At lO'* 44"^ 32.24^ 60° 23' 55" and at lO'^ 44"^ 26.17^ 60° 24' 01" (J2000), the eastern 
cloud has a mass of 7 ± 3 x 10^ Mq , a diameter of order 2 kpc, and a surface density of 
order 2 Mq pc~^. The western cloud has a mass of 4 ± 2 x 10^ Mq , a diameter <1.4 kpc, 
and a surface density of order 1 Mq pc~^. These surface densities are right at the limits 
where star formation is found to occur in other galaxies: 2-5 Mq pc~^ (Gallagher & Hunter 
1984; Skillman 1987; Taylor et al. 1994; Meurer et al. 1996; Hunter et al. 1998; van Zee et al. 
1998; Hunter et al. 2001), but much lower than would be expected for significant quantities 
of star formation (Bigiel et al. 2010). Thus, clouds like these could potentially form stars. 
Yet these clouds, at a distance of 5 disk scale lengths (9.2 kpc) from the center of the galaxy, 
are surrounded by low column density gas (1 x 10^° cm~^). 
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6. Discussion 

6.1. Radial Profiles 

Deep V-band images of 4 dwarf irregular galaxies and one BCD galaxy have been used 
to measure surface photometry down to ~ 29.5 mag arcsec"^. Deep 5-band images were also 
used for 3 of these galaxies. The surface brightness distributions in the l^-band were found 
to be exponential in radius out to 5 inner disk scale lengths (-Rd) in two cases (DDO 53, 
DDO 86), 8i?D in two other cases (NGC 4163, IZw 115), and 3Rb in one case (DDO 133). 
A clear break marks the transition from an inner exponential profile to an outer, steeper 
exponential profile in three cases. In two of these (DDO 53, DDO 86), the break occurs at 
~ 3i?D, and in the third (DDO 133) it occurs at ~ 1.8-Rd- In two cases where the B — V 
color information goes beyond the break, one (DDO 53) has a transition in the color profile 
from a slight reddening in the inner disk (from 0.12 ±0.005 to 0.36 ±0.03, a difference of 0.24 
mag) to an increasing blueness (down to —0.13 ±0.3 at the last point shown, a difference of 
0.49 mag) in the outer disk; the other (DDO 133) has a relatively flat color profile over all 
(variations up, down, up, and down again with radius covering an extreme of 0.47 ± 0.005 
to 0.30 ± 0.01, a difference of 0.17 mag). 

The remarkable regularity of these underlying exponential disks contrasts with the ir- 
regularity observed in this Hubble morphology class. Such regularity has been observed in 
the inner parts of dwarf galaxies for a long time (Patterson & Thuan 1996; Hunter et al. 
1998; Bremnes et al. 1998, 1999, 2000; Hopp 1999; van Zee 2000; Taylor et al. 2005; Hunter 
& Elmegreen 2006; Amorin et al. 2007; but for exceptions see, for example. Hunter & Gal- 
lagher 1985; Hunter & Plummer 1996; Herrmann et al., in preparation), but now we see that 
it extends into the extreme outer disk as well. Somehow the accumulated star formation 
and radial mass migration over a Hubble time converts the irregular activity at any one time 
into a standardized disk profile. This conversion seems to be done without grand design 
spiral or bar torques (only one case, DDO 133, is clearly barred and none have obvious 
spirals, but see below), and without significant shear in Hi velocity maps for at least the 
inner exponential parts. A lack of bars would mean that the Hohl (1971) mechanism for the 
formation of exponential profiles by bar and spiral torques is not operating. A lack of shear 
would suggest that the Lin & Pringle (1987) mechanism in which the star formation rate is 
proportional to the shear viscosity is also not operating. 

GALEX ultraviolet observations were also examined for these galaxies. In four cases, 
the NUV observations go out as far in the disk as the deep l^-band images. The NUV profiles 
are also exponential and have about the same scale lengths as the V-band profiles. In two of 
the three cases with exponential breaks in V-band, there is also a break in the NUV profile 
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at the same radius. In the one case with a B — V color change at the break radius (DDO 
53), there is no evident break in the NUV profile at all, even though it extends for one inner 
disk scale length beyond the V break. This lack of a NUV break in DDO 53 could be related 
to the increasing blueness with radius seen in B — V: the blueness could compensate for the 
steeper mass profile seen in the V^-band and make the NUV-band appear continuous. 

The similarity between the NUV and V^-band profiles introduces another puzzle about 
the formation and evolution of these galaxies. The average radial profile of the star formation 
rate over the last several tenths of a Gyr, as observed in the NUV, follows the underlying 
mass distribution established by star formation and radial migration over the last Gyr or 
more, as observed in V. There was no radial mixing of these disks on timescales longer than 
~ 1 Gyr because that would significantly affect the color profiles. For example, the outer 
exponential disk should be redder than the inner disk if it was formed or augmented by spiral 
arm scattering of inner disk stars ~ 1 Gyr old or more (e.g., Roskar et al. 2008). Such a 
color change has been observed for spiral galaxies (Bakos, Trujillo & Pohlen 2008), but in 
DDO 133 there is no color change and in DDO 53 the outer disk is bluer. Neither galaxy 
has spiral arms anyway. 

Consider an example using the population synthesis models of Bruzual & Chariot (2003). 
For a Salpeter IMF (the results are essentially the same for a Chabrier IMF), a burst of star 
formation older than 10^'^ yrs has NUV — V > 2 mag, and one older than 1 Gyr has 
NUV — V > 3.5 mag. Young stars have NUV — V < —1. These values were determined 
by integrating the model stellar population spectra over the filter functions used in the two 
passbands. A plot of them is shown in Figure 21 of our previous paper (Hunter et al. 2010). 
If we took a group of stars that formed 10^ years ago and moved it to a different place 
in the disk, then that group would be 3.5 magnitudes fainter in NUV than in V because 
of population fading. This 3.5 mag difference means that the displaced group may still 
stand out in V but be nearly invisible in NUV. Then the two passbands will have different 
intensity distributions. This difference of 3.5 mag is so large that even with more continuous 
star formation than a burst, and with some mixing of the group into other stars, there should 
still be a significant difference between the NUV and V intensities of stars > 1 Gyr old. Our 
observations do not show this difference, however. The similarity between the NUV and V 
profiles in our galaxies implies that there has not been much migration of stars from the 
inner disk to the outer disk on timescales longer than ~ 1 Gyr. 

The observation of similar radial profiles in NUV and V-bands (and also in J-band 
for at least the inner regions - see Hunter & Elmegreen 2006) suggests that star formation 
occurred at a steady rate. Leroy et al. (2008) also found that the star formation rate in 
dwarfs was about constant over a Hubble time. In only one of our cases (NGC 4163) is the 
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average star formation rate observed in the NUV different from the average observed in the 
l^-band by as much as a factor of 10 (the UV rate is less than the V rate at mid-disk), but in 
this case, the overall kinematics of the H I is very peculiar, suggestive of a major disruption 
in the recent past. Still, NGC 4163 has a near-perfect exponential disk in the l^-band out 
to 8 scale lengths with no break; it is the NUV profile that deviates from exponential at 
mid-radius. 

Ha observations are available for all of these galaxies from our previous surveys. The 
Ha radial profiles extend out to only 2Rd or 3i?D- In the case with the shortest l^-band 
break radius (DDO 133), the Ha extends beyond the breaks in V and NUV and also has 
a break at the same radius. In all of the galaxies, the Ha profile is roughly similar to the 
NUV profile, including local deviations from the exponential. This can be seen in Figure 



20 where we plot the ratio SFRir;7V'/SFRj:/Q,. The implication is that current star formation 
rates, seen in Ha, are about the same at each radius as the average rate over the last few 
hundred Myr. The Ha rates today are also about the same, to within a factor of ~ 3, as the 



V^-band long-term average rates in three of the galaxies (see Figure 18). The two galaxies 
where these rates differ the most are NGC 4163, which has highly disturbed kinematics as 
mentioned above (and a current rate in Ha and the UV that is much lower than the \^-band 
average at mid-radius), and the BCD galaxy, IZw 115, which also has a low Ha compared 
to V in the inner parts (although the UV rate is the same as the V rate there). IZw 115 is 
evidently post-starburst. 



6.2. Hydrogen Gas 

The final piece of evidence related to star formation and galactic structure is the gas 
distribution. For this, we examined Hi observations taken at the VLA for this project. The 
Hi radial profiles are smoothly declining but flatter than the stellar light profiles in DDO 
86, NGC 4163 and IZw 115, and about the same as the stellar light in DDO 133. In DDO 
86, there is an H I hole in the center and a ring around it that is not present in the optical. 



These trends show up in the SFR-Sgas diagrams of Figure 19 too. DDO 86, DDO 133, and 
IZw 115 have steep near-power-law relations between the \^-band or FUV SFRs and the Hi 
column densities, with a slopes of 3 ± 1 to 4 ± 1. NGC 4163 and IZw 115 have a scattering 
of points with a somewhat more organized relation for FUV than V-band SFRs. The near- 
power-laws are consistent with other studies of star formation and H I column density, which 
also show steep slopes. Bigiel et al. (2010) got a slope of 1.7 for the relation between SFR 
per unit area and H I column density in faint outer regions. Their plot of J^hi versus T^sfr 



for dwarfs in their Figure 8 is similar in shape to ours for the FUV SFR in Figure 19 Star 
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formation follows unseen CO (Bigiel et al. 2011), which is a small fraction of the Hi mass 
in these faint regions (Schruba et al. 2011). The timescale for conversion of Hi into stars is 



the ratio of the x axis values to the y-axis values in Figure 19, or approximately 10^^ years 
(see also Bigiel et al. 2010). 

The Hi distributions are clumpy with giant clouds that appear associated with star 
formation. These clouds are most obvious in the inner regions where the star formation 
rates are highest, although faint Hi clouds can also be seen in the outer parts where past- 
average star formation rates are very low. Observable clouds typically have masses of around 
10^ Mq or larger. A few faint clouds in the outer part of DDO 86 with masses of several 
xlO^ Mq were discussed as possible future sites of star formation. In general there are 
no spiral arms or systematic streaming motions in H I . The velocity fields reflect normal 
rotation curves except for NGC 4163 and DDO 53. 

The most extended Hi rotation curves are available for DDO 86, DDO 133, and IZw 
115. The first two have rising inner rotation curves that flatten out at the break radius in the 
V^-band profile. The third has a rising rotation curve in the inner part and a falling rotation 
curve in the outer part, with no break in the exponential profile for comparison. There could 
be a small companion galaxy near IZw 115 that affects the rotation curve in the outer part. 
The coincidence in two cases between the turnover point in the Hi rotation curve and the 
break point in the photometric radial profile is interesting but maybe not significant. In our 
larger sample (Hunter & Elmegreen 2006), there was no correlation between the photometric 
break radius and the rotation curve turnover radius: only one case had these features at the 
same radius, 5 had the break at a larger radius than the turnover and 4 had the break at a 
smaller radius than the turnover. In larger spiral galaxies, the exponential break radii are 
typically at ~ 2 — 5-Rd (Pohlen & Trujillo 2006) and the turnover radii are at ~ 1-Rd, so 
there is no coincidence of radii for spirals either. 

The H I rotation speeds in these dwarf galaxies are typically low: in the three cases 
with extended rotation curves, they peak at 40 km s~^ for DDO 86 and DDO 133, and at 
15 km s^^ for IZw 115. The Hi velocity dispersions are approximately constant with radius 
at around 10 km s~^. For IZw 115, the average rotation speed is comparable to the velocity 
dispersion inside SRd and possibly in the outer regions too where the rotation curve drops. 
At the peak of the rotation curve in IZw 115, the ratio of the rotation speed to dispersion is 
remarkably low, ~ 1.4. Such slow rotation gives IZw 115 very little dynamical integrity, yet 
it has a near-perfect exponential disk for 8-Rd- 

Slowly rotating galaxies should be easily perturbed by interactions with other galaxies 
and dark matter concentrations, by clumping instabilities in their disks, and by gas blowout 
during star bursts. This sensitivity to asymmetric forcing and rapidly-changing potentials 
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makes the long-range exponential structure of the disks difficult to understand, unless it 
results from such forcings. The uniform rotation observed in 2 out of 5 cases suggests that 
their recent lives have been quiescent, at least for the last several rotation periods in their 
outer regions. The rotation period at the outer radius of detection for the H I rotation curve 
in DDO 133 is 0.9 Gyr. For DDO 86, it is 1.9 Gyr. The SFR ratios shown in Figure 18 are 
roughly consistent with this long term quiescence. 



The maps of velocity dispersion, "mom2" , in Figures [13] to [17] show an inverse correla- 
tion with the maps of column density, "momO," in the sense that higher column densities 
correspond to lower velocity dispersions. We found this before for NGC 2366 (Hunter et al. 



2001). Figure 21 shows this anti-correlation better. Each panel shows a map of momO in 
red and mom2 in green. Overlapping regions turn out yellow. The maps show mostly red 
and green regions with only a little yellow. This means that the dispersion is high where the 
column density is low, and vice versa. An anticorrelation is roughly consistent with pressure 
uniformity in the sense that low density regions have higher turbulent speeds. Piontek & 
Ostriker (2005) found this in simulations of the magneto-rotational instability. 

Trujillo & Bakos (2010) looked for a correlation between possible inflection points in 
the H I velocity dispersion and exponential break radii in the disks of local spiral galaxies. 
They were considering that the breaks could result from gas disk flaring, but found nothing. 
We see nothing like that here either. 



6.3. Gravitational Stability 

The Hi appears to be gravitationally stable by itself. The ratios of the gas surface 
densities to the Toomre (1964) threshold densities are 0.6 or less in DDO 86 and DDO 
133, 0.1 or less in NGC 4163, and 0.3 in IZw 115. There are stars mixed with the gas, 
however, and the stellar velocity dispersions are probably not much larger than the gaseous 
dispersions, considering the low rotation speeds. This means that stars contribute to the 
instability in approximate proportion to their mass (Jog & Solomon 1984). 

We have estimated the mass surface densities of stars in our galaxies using the l^-band 
radial proffies and the B — V colors, where available. The calibration comes from M/L tables 
in Bell & de Jong (2000). We find that the stellar mass surface densities in the outer parts 
of our galaxies are comparable to the gas surface densities, or slightly smaller, to within a 
factor of a few. This can also be derived from the star formation rates given in Table [8] 
which are ~ lO""^'^ Mq yr~^ kpc~^ in the \^-band at fiy = 29.5 mag arcsec"^. Multiplying 
this by 10^° years and converting to pc~^, the result is a total stellar surface density today 
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of ~ 0.2 Mq pc~ in the far outer parts. Figures |8| to |12| give Sgas there, and the average is 
about 1 Mq. Thus gravitational forcing from stars may be ~ 20% of that from gas during 
the growth of an instabihty, leading to an upward correction to the surface density ratios by 
~ 20%. The correction would be larger in the inner regions, where the stellar contribution is 
larger. These corrections make DDO 86 and DDO 133 only slightly stable, which is typical 
for galaxies. For example, Leroy et al. (2008) showed that the gas+star stability parameter 
equals about 1.5 for a wide range of radii inside dwarfs and spirals. The corresponding ratio 
of column density to critical column density is the inverse of this, 0.6. This is not much 
different than it is for DDO 86 and DDO 133. 

NGC 4163 is not a good case for this stability analysis because the rotation curve is 
peculiar (§4.3). IZw 115 could also have a peculiar rotation curve because of the companion 
galaxy, but even in the inner parts where the rotation curve is rising, the H I column density 
is less than the critical value by a factor of 5. These inner parts are probably post-starburst, 
however, so they could well be very stable. 

The situation with instabilities is more complicated than this because all of these galax- 
ies have relatively thick disks, as inferred from the high ratios of HI velocity dispersion 
to rotation speed. Thickness H decreases the in-plane accelerations from self-gravity at 
wavenumber fc by a factor of ~ (1 + kH)~^ (Vandervoort 1970), which should be ~ 0.5 — 0.7 
in our case. This lowers the effective stability ratio S/S'^j.^^ by the same factor, since 
^crit — /^(^(l + kH)/7iG for Q' = 1 with thickness corrections. Here, a is the velocity disper- 
sion, S is the mass column density, and Q' = Ka{l + kH)/nGT, is the thickness-corrected 
Toomre Q parameter for epicyclic frequency n (Elmegreen 2011). The ratios derived above 
for gas+stars in DDO 86 and DDO 133 without thickness corrections are comparable to 
those given by Leroy et al. (2008) for dwarf and spiral galaxies without thickness cor- 
rections. Therefore both are likely to have comparable stability ratios when thickness is 
included, which means both are slightly more stable than without the thickness correction. 
The ratios for the outer parts of dwarfs would still be smaller than they are for spirals. 

Gravitational instabilities in dwarf irregular galaxies should look different than they 
do in large spiral galaxies. The characteristic length for spirals or clumps that form in 
a gravitational instability lies in the range a^/ (vrGE') (1 ± [1 — Q"^]) . This is centered at 
cr^/ (vrGS'), which should be compared with the disk scale length. In dwarfs with comparable 
velocity dispersions for the gas and stars, we can use the total mass column density for E' 
and include both components this way. In the outer parts of dwarf galaxies, S ~ 10 Mq 
pc~^, S' ~ 0.5S, and a ~ 10 km s~\ so a^/vrGS' = 1.4 kpc. For spiral galaxies, the length 
scale is about the same because the velocity dispersion for the dominant stars is higher by a 
factor of ~ 3, and the mass column density is higher by a factor of ~ 10. Dwarf galaxies have 
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smaller exponential scale lengths than spirals, however. Three of the dwarfs studied here 
have -Rd < 0.7 kpc and two others have Rd ~ 2 kpc (Table |4]). In spiral galaxies, Rb ranges 
between ~ 1 and 5 kpc (van der Kruit 1987), which is larger than in dwarfs by a factor of 
2 or more. This difference means that the characteristic length of a significantly gravitating 
structure in a large spiral galaxy readily fits into the disk, allowing coherent structures like 
spiral waves to form. Most dwarfs are too small for this. In addition, the bright parts of 
dwarf galaxy disks occur in the rising parts of their rotation curves, but the main parts of 
spiral galaxy disks are in the fiat parts of their rotation curves. Gravitational perturbations 
therefore tend to grow as unsheared clumps in the dwarfs, but they produce spiral wakes 
and long sheared arms in spirals. 

Some dwarf galaxies have spirals, but this situation is rare (Schombert et al. 1995). A 
dwarf elliptical in the Virgo cluster, IC 3328, has a weak spiral, but it may be confined to a 
thin sub-component of the main stellar disk where the gravitational length is short (Jerjen 
et al. 2000; De Rijcke et al. 2003; Lisker & Fuchs 2009). Other dwarf spirals in Virgo were 
observed by Barazza et al. (2002). 

The same considerations apply to the far-outer parts of spiral galaxies, where S is 
comparable to that in dwarfs. Then the characteristic length for instabilities may be 10 
times higher than it is in the inner spiral disk, and locally much larger than the scale length. 
Thus the outer parts of spiral disks should have difficulty generating self-gravitating waves 
on their own, like dwarf galaxies. However, the outer disks of spiral galaxies can have spiral 
waves that propagate there from the inner disk (Bertin & Amorisco 2010) or were resonantly 
excited by perturbations from the inner disk (Debattista et al. 2006). Dwarf galaxies have 
neither source of waves and should produce only clumps of self-gravitating gas. 

Disk stability is also relevant to the origin of double exponential profiles, according to 
the model by Debattista et al. (2006). This model involves angular momentum transfer 
around bars and spirals. Debattista et al. (2006) showed that the inner parts of barred 
galaxies with single exponential profiles get fiatter with time, out to the outer Lindblad 
resonance of the spiral that is driven by the bar. The original exponential disk then shows 
up as the outer steep part of the double exponential, and the bar-modified part is the inner 
fiat region. The degree of fiattening depended on Q in these models. Disks that were more 
gravitational stable, with higher Q, did not evolve as much as less stable disks. This result 
suggests that dwarf galaxies that are weakly self-gravitating may not be responsive to bar 
torques that modify a single exponential into a double exponential. 

Our galaxies with the highest ratios of Sgas/Scrit, namely DDO 86 and DDO 133, also 
have the most obvious breaks in their exponential profiles. DDO 133 is a barred galaxy 
too. However, the break radius in DDO 133 is at the end of the bar and not the end of a 
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spiral as in the Debattista et al. model. Also in DDO 86, the star formation region and 
inner bright core in Figure |2] is elongated in the east-west direction, reminiscent of a bar. 
The end of this elongation corresponds to the break radius. Similarly in DDO 53, which is 
the only other case here with a double exponential, there is an elongation of the main disk 
star-forming region along a position angle of ~ 150°, which is transverse to the major axis. 
The deprojected length of this elongation is comparable to the break radius. There is not 
enough H I data to determine Egas/Ecrit for this galaxy. These observations suggest that the 
inner exponentials of our galaxies were made flatter than the outer exponentials by mass 
transfer connected with a bar. There are no spirals here, so the bar may be doing here what 
the spirals did in Debattista et al. (2006). 

We looked for a correlation between photometric breaks and the presence of bars in our 
larger sample (Hunter & Elmegreen 2006). Out of 22 dwarf galaxies with double exponentials 
of the type discussed here, 8 have some evidence for a bar (even if they are not classified as 
Hubble bar types), and the rest do not. DDO 133, studied here, was one of the barred cases 
with a double exponential recognized in that previous study. Among the barred dwarfs, the 
ratio of the break radius to the end of bar radius is 1.8 ± 0.6 on average. Since there are 
breaks without bars and bars without breaks, we do not see a clear correlation between the 
two at this time. 



6.4. Two Peculiarities 

We summarize this discussion by highlighting two peculiar things about the present 
observations: the extraordinarily regular exponential disks in dwarf irregular galaxies, and 
the continuity of star formation throughout these disks and over time. Exponential disks are 
made in part by cosmological infall (Freeman 1970; Fall & Efstathiou 1980; van der Kruit 
1987) and subsequent dynamical torques from self-gravitating disk structures. There is no 
regularity to these structures in dwarfs and no systematic presence of bars or spirals. There 
are only clumps of star formation and maybe some old disk stars mixed in with those clumps. 
It is conceivable that the clumps of star formation make or maintain the exponential disk, 
as in numerical simulations of equally clumpy galaxies in the early universe (Bournaud et 
al. 2007). Or, there could be bar-like dark matter halos (Bekki & Freeman 2002; Tutukov 
& Fedorova 2006) that redistribute angular momentum, while also driving in gas to fuel 
occasional low- level bursts (Hunter & Elmegreen 2004). 

As for star formation itself, this complex process involving mostly gas seems to follow the 
underlying exponential disk of old stars better than one would expect from such a stochastic 
process. Leroy et al. (2008) found the same thing for dwarfs and the outer parts of spiral 
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galaxies. Molecular cloud distributions follow the stellar exponential even when the molecules 
are an extremely low fraction of the gas mass (Schruba et al. 2011). How do the molecules 
know to do this? In the current paradigm, star formation follows only the molecules, and 
molecular abundances depend mostly on the pressure (e.g., Bigiel et al. 2008), but why does 
this end up reinforcing the underlying exponential disk? There is a purely radial dependence 
to the cloud and star formation rates, in addition to a gas dependence, as in the dynamical 
law of star formation, where the rate is determined by the orbit time (Kennicutt 1998). Bigiel 
et al. (2010, Figure 7) also found that for a given H I column density, the star formation 
rate decreases with increasing radius in spiral and dwarf galaxies. The problem with this is 
that orbital timescales usually vary with radius as a power law, not an exponential. 

Ostriker et al. (2010) proposed that the star formation rate scales with the total gas 
surface density multiplied by the square root of the midplane density from stars and dark 
matter. They get good agreement between this prediction and the radial profiles of star 
formation in several spiral galaxies. The present results suggest that if this quantity applies 
to dwarf irregulars, then it should also be proportional to the column density of existing 
stars so that the color profile does not change much over at least a Gyr. In fact they find the 
star formation rate per unit gas mass proportional to the surface density of existing stars 
and dark matter in the case where the vertical velocity dispersion is constant (because then 
the midplane density scales with the square of the stellar and dark matter column density). 
This matches our result fairly well because the gas column density varies weakly with radius. 

Another model is that in highly stable regions, molecular cloud formation and star 
formation is somehow triggered by field stars. Perhaps Type la supernovae are important, 
or the turbulence that is driven by these supernovae and other field-star pressures. Star 
formation would then follow the underlying stellar density as long as there is enough gas 
to compress. Since the rate is approximately constant over time, sequential triggering by 
young star pressures could also be important. Then star formation would move around as 
illustrated by Gerola et al. (1980). A way to check these models would be to find triggering 
sources near young stars. Also, the star formation rate should be more tightly correlated 
with the stellar column density than with the gas column density, as long as there is enough 
gas to be triggered. 

The break radii in exponential disks could be a clue to these processes if the properties of 
the gas and cosmic rays near regions of star formation were found to differ inside and outside 
of this radius. In the present sample of galaxies, we found no clear correlation between the 
existence of a break and other properties, such as bars, spirals, color jumps, or the turnover 
in the rotation curve. We found no such correlations in our bigger sample either (Hunter & 
Elmegreen 2006). 
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7. Summary 

We present ultra-deep V, B, and GALEX UV images obtained in order to detect and 
measure the extreme outer stellar disks of a sample of 4 dim galaxies and one BCD galaxy. 
Our l^-band surface photometry extends to 29.5 magnitudes arcsec"^. We convert our new 
FUV and V-band photometry, along with Ha photometry obtained in a larger survey into 
radial star formation rate profiles that are sensitive to different timescales — from 10 Myrs 
to the lifetime of the galaxy. We also present H I -line emission data and compare the stellar 
distributions, surface brightness profiles, and star formation rate profiles to H I -line emission 
maps, gas surface density profiles, and gas kinematics. Our data lead us to two remarkable 
general observations. 

First, the exponential disk in these dwarf irregular galaxies are extraordinarily regular. 
We observe that the stellar disks continue to decline exponentially as far as our measurements 
extend. In addition, as is common in outer spiral and dwarf disks, the slope of the l^-band 
surface brightness profiles change abruptly in the outer disks in most of our systems. In three 
of our galaxies the outer profile becomes steeper, but the break in the \^-band surface bright- 
ness profile is mimicked in the NUV profile for two of these systems. The galaxies are lumpy 
with star formation and the Hi distributions and kinematics have significant non-ordered 
motions. Yet, the regularity of the exponential profiles in spite of these peculiarities indi- 
cates that the sporadic processes that have built the disks — star formation, radial movement 
of stars, and perhaps even perturbations from the outside — have, nevertheless, conspired to 
produce standard disk profiles. 

Second, there is a remarkable continuity of star formation throughout these disks over 
time. In four out of five of our galaxies the star formation rate in the outer disk measured 
from the FUV tracks that determined from the V"-band, to within factors of a few (<5), 
although the FUV-derived rate refiects a shorter timescale than the V-band derived rate. 
This requires star formation at a fairly steady rate over the galaxy's lifetime. But how the 
galaxy achieves this is not clear. The Hi surface density profiles generally decline with 
radius more shallowly than the stellar light and, hence, star formation rate profiles. The star 
formation instead appears to follow the older stars. There is a general correlation between 
azimuthally-averaged log T^gas and log SFR, similar to what is seen by others, except for some 
radii in the two galaxies with the disturbed kinematics. The gas is marginally gravitationally 
stable against collapse into clouds in these galaxies. 

Outer stellar disks are being traced to ever fainter levels in dwarf galaxies, and these 
observations are challenging our concepts of star formation and disk growth over cosmic 
times. Clearly, understanding the processes that mold outer disks, as well as understanding 
how stars form in the extreme environment of outer disks, are important to our understanding 
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of galaxy formation and evolution. 
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Fig. 1. — Left and Upper Right: V, Ha, and NUV images of DDO 53, shown as the logarithm 
in order to show a wider range of brightness levels. The white ellipses show the breaks in 
the surface photometry and the outer-most ellipse that was used in the surface photometry. 
Bottom Right: l^-band image with contours of integrated Hi emission superposed: 1.9 to 



64.3 in steps of 5.67x10 cm . The large red ellipse denotes the radius of the major break 
in the ^-band surface photometry. The small red ellipse in the lower left corner shows the 
FWHM of the VLA beam. 
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Fig. 2. — As in Figure 1, for DDO 86. The Ha image is not shown as the logarithm. For 
DDO 86, photometry was measured at the position angle defined by the moment zero map 
fit to the Hi kinematics. Contours of integrated Hi emission are 13.5 to 22.9 in steps of 



5.38x1020 cm-2. 



41 



DDO 133 - 


LOG U 




• 




* 

E 


/ 1 ■ 




• 


^ 


« 


*lJI 


i * h 


»+ 






tfltfg/ii 


M )i 






■ . « 


JbHb 


1^ 




• 


* < ' 


. \-' • ■ . 








r 


\, • ■ 


* 


* ' ■ 




■ . 50" .' 


■"----^ - 


<«H 


4 




* 
* • 


• 




* 


* 





Fig. 3. — As in Figure 1, for DDO 133. The Ha image is not shown as the logarithm. 
Contours of integrated Hi emission are 0.06, 0.3, 1.58 to 11.7 in steps of 1.26x10^ 



120 ^^-2_ 
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Fig. 4. — As in Figure 1, for NGC 4163. The Ha image is not shown as the logarithm. 
Contours of integrated Hi emission are 1.8 and 3.03 to 7.89 in steps of 2.43x10^^ cm~^ and 
1.7 to 8.1 in steps of 0.91x10^° cm'^. 
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Fig. 5. — As in Figure 1, for IZw 115. The Ha image is not shown as the logarithm. 



Contours of integrated Hi emission are 1.3 to 27.2 in steps of 2.59x10 cm~ . The region 
of stars that is detached from the main body of IZw 115 is marked in the V and NUV images. 
The southern H I tail in the color composite image extends over the detached region. 
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Fig. 6. — Deep V, and in some cases B — V, surface photometry with different methods 
of determining the background/foreground object contamination: "edit" refers to editing 
and masking individual unwanted objects, "stat" refers to using the surrounding field to 
determine a statistical correction, and "stat-2" refers to the statistical method carried out 
by a different co-author. We also compare the deep surface photometry to the shallower 
survey surface photometry from Hunter & Elmegreen (2006). Uncertainties are only shown 
where they are larger than the symbols. R25 and Rh mark the radius at the i?-band 25 mag 
arcsec"^ isophote and the Holmberg radius {B =26.66 mag arcsec"^), respectively. 
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Fig. 7. — Deep GALEX NUV and FUV — NUV surface photometry with different meth- 
ods of determining the background/foreground contamination: "edit" refers to editing and 
masking unwanted objects, "stat" refers to using the surrounding field to determine a sta- 
tistical correction, and "psf ' refers to point spread function (PSF) fitting and subtraction. 
Uncertainties are only shown where they are larger than the symbols. i?25 and Rh mark the 
radius at the 5-band 25 mag arcsec"^ isophote and the Holmberg radius {B =26.66 mag 
arcsec"^), respectively. 
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Fig. 8. — Azimuthally- averaged surface photometry and colors of DDO 53, corrected for 
reddening, plotted against radius from the galaxy center normalized to the disk scale-length 
Rd- Top: Hi properties, including the gas surface density, rotational velocity, velocity 
dispersion of the rotational component, and the ratio of gas surface density to critical density. 
Tjcrit is the critical threshold density for large-scale gravitational instabilities in a differentially 
rotating thin disk (Toomre 1964). S^as includes H I and He. Bottom: Optical and ultraviolet 
properties, including V, Ha, and NUV surface photometry. i?25 and Rh mark the radius 
at the 5-band 25 mag arcsec"^ isophote and the Holmberg radius (26.66 mag arcsec"^), 
respectively. The dashed vertical black lines mark breaks in the \^-band profile. 
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Fig. 9. — Azimuthally- averaged surface photometry and colors of DDO 86, corrected for 
reddening, plotted against radius from the galaxy center normalized to the disk scale-length 
Rd- Top: Hi properties, including the gas surface density, rotational velocity, velocity 
dispersion of the rotational component, and the ratio of gas surface density to critical density. 
Ticrit is the critical threshold density for large-scale gravitational instabilities in a differentially 
rotating thin disk (Toomre 1964). S^as includes H I and He. Bottom: Optical and ultraviolet 
properties, including V ^ Ha, and NUV surface photometry. i?25 and Rh mark the radius 
at the 5-band 25 mag arcsec"^ isophote and the Holmberg radius (26.66 mag arcsec"^), 
respectively. The dashed vertical black line marks the break in the \^-band profile. 
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Fig. 10. — Azimuthally-averaged surface photometry and colors of DDO 133, corrected 
for reddening, plotted against radius from the galaxy center normalized to the disk scale- 
length Rd. Top: H I properties, including the gas surface density, rotational velocity, velocity 
dispersion of the rotational component, and the ratio of gas surface density to critical density. 
Ticrit is the critical threshold density for large-scale gravitational instabilities in a differentially 
rotating thin disk (Toomre 1964). S^as includes H I and He. Bottom: Optical and ultraviolet 
properties, including V ^ Ha, and NUV surface photometry. i?25 and Rh mark the radius 
at the 5-band 25 mag arcsec"^ isophote and the Holmberg radius (26.66 mag arcsec"^), 
respectively. The dashed vertical black line marks the break in the \^-band profile. 
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Fig. 11. — Azimuthally-averaged surface photometry and colors of NGC 4163, corrected 
for reddening, plotted against radius from the galaxy center normalized to the disk scale- 
length Rd. Top: H I properties, including the gas surface density, rotational velocity, velocity 
dispersion of the rotational component, and the ratio of gas surface density to critical density. 
Ticrit is the critical threshold density for large-scale gravitational instabilities in a differentially 
rotating thin disk (Toomre 1964). S^as includes H I and He. Bottom: Optical and ultraviolet 
properties, including V ^ Ha, and NUV surface photometry. i?25 and Rh mark the radius 
at the 5-band 25 mag arcsec"^ isophote and the Holmberg radius (26.66 mag arcsec"^), 
respectively. The dashed vertical black line marks the break in the \^-band profile. 
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Fig. 12. — Azimuthally-averaged surface photometry and colors of IZwllS, corrected for 
reddening, plotted against radius from the galaxy center normalized to the disk scale-length 
Rd- Top: Hi properties, including the gas surface density, rotational velocity, velocity 
dispersion of the rotational component, and the ratio of gas surface density to critical density. 
Ticrit is the critical threshold density for large-scale gravitational instabilities in a differentially 
rotating thin disk (Toomre 1964). S^as includes H I and He. Bottom: Optical and ultraviolet 
properties, including V ^ Ha, and NUV surface photometry. i?25 and Rh mark the radius 
at the 5-band 25 mag arcsec"^ isophote and the Holmberg radius (26.66 mag arcsec"^), 
respectively. The dashed vertical black line marks the break in the \^-band profile. 
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Fig. 13. — Moment maps and a major-axis position-velocity diagram of DDO 53. (a) 
Integrated Hi map (moment 0). Contours run from to 2 Jy beam^^ km s^^ with a spacing 
of 0.2 Jy beam~^ km s~^. The RA,DEC axes can be used to match the field of view with 
that in Figure [I] {middle right panel), (b) Intensity- weighted mean velocity field (moment 
1). Contours run from —20 km s~^ to 35 km s~^ with a spacing of 2 km s~^. (c) Hi velocity 
dispersion map (moment 2). Velocity contours run from 1 to 17 km s~^with a spacing of 
2 km s~^. (d) Hermite h^ velocity field. Contours run from —20 km s^^to 35 km s~^ with 
a spacing of 2 km s~^. (e) Position-velocity diagram taken along the average position angle 
(155°) of the major axis. Contours start at +3cr in steps of 3cr. The dashed lines indicate the 
systemic velocity and position of the kinematic center derived from the tilted-ring analysis. 
Overplotted is the hermite velocity field rotation curve corrected for the average inclination 
(27°). 
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Fig. 14. — Moment maps and a major-axis position-velocity diagram of DDO 86. (a) 
Integrated Hi map (moment 0). Contours run from to 0.38 Jy beam~^ km s~^ with a spacing 
of 0.08 Jybeam"^km s~^. The RA,DEC axes can be used to match the field of view with 
that in Figure |2] {middle right panel), (b) Intensity- weighted mean velocity field (moment 
1). Contours run from 1000km s^^ to 1050km s~^ with a spacing of 5km s^^. (c) Hi velocity 
dispersion map (moment 2). Velocity contours run from 1 to 20 km s~^with a spacing of 
2 km s~^. (d) Hermite h^ velocity field. Contours run from 1000 km s^^ to 1050 km s~^ with 
a spacing of 5 km s~^. (e) Position-velocity diagram taken along the average position angle 
(317°) of the major axis. Contours start at +3a in steps of 3a. The dashed lines indicate the 
systemic velocity and position of the kinematic center derived from the tilted-ring analysis. 
Overplotted is the hermite velocity field rotation curve corrected for the average inclination 
(32°). 
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Fig. 15. — Moment maps and a major-axis position-velocity diagram of DDO 133. (a) 
Integrated Hi map (moment 0). Contours run from to 3 Jybeam~^km s~^ with a spacing 
of 0.5 Jybeam~^km s~^. The RA,DEC axes can be used to match the field of view with 
that in Figure |3] {middle right panel), (b) Intensity- weighted mean velocity field (moment 
1). Contours run from 280km s~^to 400 km s~^ with a spacing of 5 km s~^. (c) Hi velocity 
dispersion map (moment 2). Velocity contours run from 1 to 20 km s~^ with a spacing of 
2 km s~^. (d) Hermite h^ velocity field. Contours run from 280 km s~^ to 400 km s~^with 
a spacing of 5 km s~^. (e) Position-velocity diagram taken along the average position angle 
(358°) of the major axis. Contours start at +lcr in steps of 3cr. The dashed lines indicate the 
systemic velocity and position of the kinematic center derived from the tilted-ring analysis. 
Overplotted is the hermite velocity field rotation curve corrected for the average inclination 
(46°). 



54 



NGC 4163 




36° 7 



a oh^ rjm^ c-S -|imcrc:S +0,04 +0.02 -0.02 -0,04 

Offset (orcmin) 



RA (2000.0) 



Fig. 16. — Moment maps and a major-axis position- velocity diagram of NGC 4163. (a) 
Integrated Hi map (moment 0). Contours run from to 1.9 Jybeam~^km s^^ with a 
spacing of 0.2 Jybeam~^km s~^. The RA,DEC axes can be used to match the field of 
view with that in Figure |4] {middle right panel), (b) Intensity-weighted mean velocity field 
(moment 1). Contours run from 150km s~^to 170km s~^with a spacing of 2km s~^. (c) 
Hi velocity dispersion map (moment 2). Velocity contours run from 1 to 17km s~^with a 
spacing of 2 km s~^. (d) Bulk velocity field. Contours run from 150 km s~^ to 170 km s~^ with 
a spacing of 2 km s~^. (e) Position-velocity diagram taken along the average position angle 
(15°) of the major axis. Contours start at +la in steps of 3a. The dashed lines indicate the 
systemic velocity and position of the kinematic center derived from the tilted-ring analysis. 
Overplotted is the hermite velocity field rotation curve corrected for the average inclination 
(48°). 
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Fig. 17. — Moment maps and a major-axis position- velocity diagram of IZw 115. (a) Inte- 
grated Hi map (moment 0). Contours run from to 0.2 Jybeam~^km s^^ with a spacing 
of 0.03 Jybeam~^km s~^. The RA,DEC axes can be used to match the field of view with 
that in Figure |5] {middle right panel), (b) Intensity- weighted mean velocity field (moment 
1). Contours run from 600km s~^to 700 km s~^ with a spacing of 5 km s~^. (c) Hi velocity 
dispersion map (moment 2). Velocity contours run from 1 to 15 km s~^with a spacing of 
2 km s~^. (d) Hermite h^ velocity field. Contours run from 600 km s^^to 700km s~^with 
a spacing of 5 km s~^. (e) Position-velocity diagram taken along the average position angle 
(51°) of the major axis. Contours start at +3a in steps of 3a. The dashed lines indicate the 
systemic velocity and position of the kinematic center derived from the tilted-ring analysis. 
Overplotted is the hermite velocity field rotation curve corrected for the average inclination 
(61°). 
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Fig. 18. — Azimuthally-averaged star formation rates (SFRs) determined from the FUV 
luminosity SFRpuv {blue) and from the Ha SFRua {red) divided by SFRy for the 5 
galaxies in our sample. SFRy is determined from the mass in starscalculated from My 
and a stellar M/Ly ratio that depends on {B — V)q, and assumes a constant SFR over 12 
Gyr. SFRfja is determined from the Ha luminosity. SFRho tells us of the past 10 Myrs; 
SFRpuv integrates over a longer time scale; and SFRy is most sensitive to the past 1 Gyr 
for constant SFRs. The horizontal dashed lines mark equal SFRs, or ratios of 1. The ratios 
are plotted as a function of the radius normalized to the disk scale length measured from 
V^-band surface brightness profiles R^. 



57 



u 

Q. 



o 

2 



a. 
to 



o 



-2 - 
-3 
-4 - 



- Ha 



r 

FUV 



-D53 
086 

-D133 
N4163 
IZw115 



- V 




Fig. 19. — Azimuthally-averaged logEga^ versus log SFR for three star formation measures: 
SFRiT'iyv, SFRhq,, and SFRy. The sohd, slanted hues are hues of gas depletion times; the 
lines are labeled with the gas depletion timescale in years. 
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Fig. 20. — Ratio of azimuthally-averaged SFR determined from the FUV luminosity, 
SFRpuv, to that determined from the Ha, SFRh^- SFRua tells us of the past 10 Myrs; 
SFRpuv integrates over several hundred million years. The horizontal dashed lines mark 
equal SFRs, or ratios of 1. The ratios are plotted as a function of the radius normalized to 
the disk scale length measured from l^-band surface brightness profiles R^. 
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Fig. 21. — Integrated Hi maps ("inoinO") (red) superposed with velocity dispersion maps 
("mom2") (green). 
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Table 5. GALEX Observations 



FUV Exp (s) NUV Exp (s) Tile Name 



DDO 53 


9419 


8595 


DDO 86 


14050 


15194 


DDO 133 


9841 


18287 


IZwllS 


8234 


7616 


NGC 4163 


10173 


12392 



GI4_015001_DDO53 

LOCK_20 

GI4_015003_DDO133 

GI4_015005JZW115 

GI4_015004_NGC4163 
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Table 6. VLA Observations 









Time on-source 


Beam size 


Galaxy 


Date 


Array Config 


(hours) 


(arcsec) 


DDO 53 


1997 Sep & Nov 


CS,D 


3.7,3.6 


28.5" X 26.7" 


DDO 86 


2006 Dec 


cs 


4.7 


15.7"xl3.5" 


DDO 133 


2004 Jun & Aug 


D 


8.2 


48.9" X 45.9" 


NGC 4163 


2004 Jun & Aug 


D 


7.5 


50.8" X 46.0" 


IZwll5 


2008 Mar 


CS 


3.4 


15.1"xl4.6" 
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